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Abstract. We summarize the x-ray properties of the complete samples 
of millisecond pulsars (MSPs) detected in our Chandra observations of 
the globular clusters 47 Tuc and NGC 6397. The 47 Tuc MSPs are pre¬ 
dominantly soft sources suggestive of thermal emission from the neutron 
star polar cap and have x-ray luminosities in a surprisingly narrow range 
(La; ~l-4 xl0^° erg s“^). The single MSP in NGC 6397 is both hard and 
apparently extended, probably due to shocked hot gas evaporating from 
its main sequence companion. In contrast to MSPs in the field and the 
cluster M28, which show correlation between x-ray luminosity and spin- 
down luminosity Lx oc E ^ with /3 ~1 - 1.4, the 47 Tuc (and NGC 6397) 
sample display a relatively tight correlation with (3 = 0.5 ± 0.15. The 
correlations of La, vs. P/2P and light cylinder magnetic field values are 
also different. It is possible the magnetic field configuration has been 
altered (by episodic accretion) for old MSPs in dense cluster cores. 


1. Introduction 

X-ray studies of millisecond pulsars (MSPs) can constrain fundamental prop¬ 
erties of their emission regions and, when combined with radio timing studies, 
their underlying neutron stars (NSs). In globular clusters both MSPs and low 
mass x-ray binaries (LMXBs), their likely progenitors, are significantly enhanced 
(per unit mass) over their values in the galactic disk by stellar and binary inter¬ 
actions. The dense cluster (core) environment needed for their excess formation 
may also alter their evolution. Thus cluster vs. field MSPs, as studied in x-rays 
and radio, can constrain intrinsic vs. extrinsic (evolutionary) properties of these 
oldest NS systems. 

We have conducted a deep Chandra survey for MSPs as well as quiescent 
LMXBs and cataclysmic variables (GVs) in the globular clusters 47 Tuc (Grind- 
lay et al. 2001a; GHEOla) and NGC 6397 (Grindlay et al. 2001b; GHEOlb). 
The full details of the MSP survey are given in Grindlay et al. (2001c; GCHOl). 
Here we present the highlights of this study, focusing on just the x-ray proper¬ 
ties of the 16 MSPs with radio timing positions in 47 Tuc (Freire et al. 2001a, 
Freire 2001) and the one in NGC 6397 (D’Amico et al. 2001; DPM) as well as 
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their comparison with the field MSP population (cf. Becker & Trumper 1997, 
1999; BT97, BT99). We defer to the full paper the discussion of the total MSP 
populations and spatial distributions, which probe cluster dynamics. 


2. X-ray Colors and Emission Models 

The 47 Tuc MSPs were found initially (GHEOla) to be soft sources. In GCHOl 
we give the detected counts in 3 bands: softcts (0.2-1 keV), mediumcts (1-2 keV) 
and hardcts (2-8 keV) for each of the 14 resolved MSPs, with counts for 47 Tuc- 
G and -I (unresolved) estimated. From these bands, we form the hardness ratios 
HRl = mediumcts/softcts and HR2 = hardcts/mediumcts and plot the MSPs, 
with counting statistics errors, in the color-color diagram shown in Figure 1 
(left). The MSP colors are clustered in a relatively narrow range of HRl and 
HR2 with 47 Tuc-J clearly harder, as was evident in the Xcolor distributions in 
GHFOla. 

Using the PIMMS tool, we construct values of HRl and HR2 for 3 simple 
models: thermal bremsstrahlung (TB), blackbody (BB) and power law (PL), 
with index values (kT or photon index) given in the caption of Figure 1 (left). 
The observed range of HR1-HR2 is roughly consistent with TB spectra with kT 
~lkeV, BB spectra with kT ~0.2-0.3keV (except for 47 Tuc-J) or PL spectra 
with photon index ~3. The weighted mean colors for all but 47 Tuc-J are 
consistent with a BB spectrum with kT ~0.22 keV, giving x-ray luminosities 
Ijx (0.5-2.5keV) ~l-4 xlO^'^ erg s~^and thus mean bolometric ^ix—boi — 2.6 x 1030 
erg s“^. 

The x-ray colors rule out TB models (surrounding column densities incon¬ 
sistent with the MSP dispersion measures; DM) and PL fits (spectral indices im¬ 
plausible). Simple BB fits for hx-boi give emission radii of only ~0.1km whereas 
H (or He)-atmosphere models (Rajagopal & Romani 1996) typically give temper¬ 
atures reduced (from BB) by a factor of ~2 and thus radii increased to ~0.4km. 
Either case suggests soft x-ray emission from a region smaller than the entire 
polar cap, as predicted in recent models of Harding & Muslimov (2001) for polar 
cap heating. Although the 3.2s temporal resolution of Chandra-AClS prevents a 
pulsation analysis, the small thermal emission area suggests the emission would 
be pulsed, with a sinusoidal pulse shape appropriate to the fractional visibility 
of the isotropically radiating thermal polar cap. In contrast, the narrower pulse 
duty cycles of ~10% for some field MSPs (and one in the globular cluster M28; 
BT99) are probably due to non-thermal beamed emission. 


3. X-ray vs. Pulsar Spin Properties 

A key question for this rich Chandra dataset is the correlation of x-ray luminosity 
and pulsar spindown luminosity E, which is found for field MSPs (with much 
more uncertain distances) to scale as Lx (0.1-2.4keV) 10 ^ E (BT97) and with 

a possibly steeper logarithmic slope (1.4) for Lx in the 2-lOkeV band (Possenti 
et al. 2001; PGG). We derive instrinsic period derivatives, Pi, corrected for 
the cluster acceleration by estimating the 3D positions of each MSP in the 
cluster from the observed DM value and the observed hot gas and thus electron 
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Figure 1. Left: X-ray color-color diagram for MSPs in 47 Tuc (solid 
diamonds) vs. tracks for models with parameters changing from lower 
left to upper right; blackbody (BB; open squares), with kT = 0.2, 0.25, 
0.3keV; power law (PL; open triangles), with photon index a = 3, 2, 
1.5, 1; and thermal bremsstrahlung (TB; open circles), with kT = 1, 
2, 3, 6keV. Weighted mean colors for the MSPs in 47 Tuc with and 
without 47 Tuc-J shown as the solid square and triangle, respectively. 
Right: Lx vs. E for MSPs in 47 Tuc (diamonds) and field (triangles). 
MSP in NGC 6397 (box) is labeled as is J in 47 Tuc. 


density in the cluster (Freire et al. 2001b) and then subtracting the cluster 
acceleration using a King model with cluster parameters derived by Meylan & 
Mayor (1986). Using a standard NS moment of inertia I = 10‘^^gcm^, we then 
derive E = for each MSP and plot them vs. Lx (0.5-2.5keV) in 

Figure 1 (right). Uncertainties in the E values are typically 0.2-0.5 in the log 
but are not shown for clarity; uncertainties in log(La,) are typically 0.2, and 
extrapolating to the ROSAT band, 0.1-2.4keV, would increase log(La;) only by 
0.1. For comparison with 47 Tuc, we plot the MSP in NGC 6397 (GHEOlb), for 
which the E uncertainty is small, and updated values (cf. GCHOl) for the 10 
field MSPs previously detected in x-rays as well as in the globular cluster M28. 

Whereas the MSPs in the field and M28 show (Figure 1, right) a correlation 
logLj, (0.1-2.4keV) = (1.13±0.15)logP - 7.5±5, the MSPs in 47 Tuc appear to 
have a weaker dependence: logLa, (0.5-2.5keV) = (0.47 ± 0.10)logP -|- 14.2 ± 3.4 
for the nominal cluster model with central velocity dispersion (T„o= 11.6 km s“^, 
where the errors (Tlo") in both correlations are due to just the scatter in the 
points. Allowing for uncertainties in the cluster model and distance gives slope 
0.48 ± 0.21 and intercept 13.8 ± 7.5. Including the errors for the E values 
estimated for the 47 Tuc MSPs, but with the approximation that unequal errors 
(on E) are simply averaged (which biases the slope to steeper values, since the 
unequal errors are much larger for smaller values of E), increases the logarithmic 
slope to (3 = 0.62 ±0.29 and offset to 9.0 ± 10.8. The best (median) estimate for 
the 47 Tuc MSPs alone is thus (3 ~0.55±0.2. Apart from the uncertain detection 
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Figure 2. Left: Lx vs. characteristic age {P/2Pi) for MSPs in 47 Tuc 
(diamonds) and field (triangles). MSP in NGC 6397 (box) is labeled 
as is J in 47 Tuc. Right: Lx vs. Bic, inferred magnetic field strength at 
the light cylinder (ric = cPI2tt) for MSPs in 47 Tuc (diamonds) and 
field (triangles). MSP in NGC 6397 (box) is labeled as is J in 47 Tuc. 
Values for MSPs in field and M28 from data given in GCHOl, with B\c 
derived assuming standard 1/r^ dependence of dipole field (see text). 
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of 47 Tuc-C (GCHOl), the MSPs in 47 Tuc have Lx (0.5-2.5 keV) values within 
a factor of ~4 despite a range of ~25 in E. Figure 1 (right) shows that 6397-A, 
the MSP in NGC 6397, is consistent with the Lx - E correlation shown by the 
47 Tuc MSPs. Including 6397-A in the fit yields f3 = 0.50 ± 0.08 (scatter only) 
or 0.55±0.15 (with averaged Tier errors for the 47 Tuc E values). 

The smaller /3 values for 47 Tuc and NGC 6397 may be due to the different 
formation histories and possibly different physical parameters of their MSPs vs. 
the field objects. In Figure 2 (left) we plot Lx vs. spindown ages, Pf2Pi, for 
MSPs in the field (and M28) vs. the 47 Tuc and NGC 6397 sample. Error bars 
on the age parameter are not shown, for clarity, but are typically ~lo'i in the log 
and primarily due to the 47 Tuc acceleration model. Despite the uncertainties, 
the correlation is striking: the field MSPs show a declining Lx with “age,” and 
the 47 Tuc MSPs appear to fall on this trend but with a flatter Lx vs. age 
slope. However, spindown ages correspond only approximately to actual ages, 
and while these estimates are consistent with the formation of most 47 Tuc MSPs 
early in the cluster history, this would not by itself provide an explanation for 
the different Lx vs. E correlations found for MSPs in 47 Tuc and NGC 6397 vs. 
the field. 

Another possible physical difference between the 47 Tuc and field MSPs 
might be magnetic field strength at the light cylinder (at which the corota¬ 
tion speed equals c), since this likely affects the relative importance of non- 
thermal (magnetospheric) emission. For an assumed dipole field, this is given 

by Bic = 9.35 X 10^Pmsec G, with P* in units of 10“^'^. In Figure 2 (right) 
we plot Lx vs. Pic. Again, considered as a homogeneous group, the MSPs in 
the field (and M28) lie on a steeper slope than the 47 Tuc MSPs. In support 
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of the hypothesis that B\c, rather than the field -Bsurf = 3.2 x 10^® G 

at the NS surface, is more closely related to Lx is the fact that the correlation 
of hx with Psurf is less defined and with larger scatter, and with logarithmic 
slopes differing even more: 0.05±0.27 for the 47 Tuc MSPs vs. 2.80±0.99 for 
the field MSPs. We note that 3 out of the 4 MSPs with Pic ^ 10^'^ G display 
x-ray emission that is seemingly magnetospheric, with the nature of emission 
from the 4th (the eclipsing MSP B1957-I-20) indeterminate (BT99; Takahashi et 
al. 2001). Conversely, field pulsars with Pic < 10^ G have x-ray emission that is 
typically either thermal or of indeterminate character (BT99). Considering the 
small numbers of such pulsars studied, and that most of them have highly un¬ 
certain distances, it seems possible that field pulsars with Pic ^ 10^ G may show 
an La;-Pic trend that is fairly flat and roughly consistent with the better deter¬ 
mined relation for the 47 Tuc MSPs. However for this interpretation to hold, a 
few field MSPs with relatively well determined distances (e.g. 0437—4715 and 
1744—1134; cf. GCHOl) must be accounted for and the even larger deviations of 
0751-1-1807 and 1024—0719 from the 47 Tuc correlation line would require factor 
^3 adjustments in these MSP distances. 


4. Conclusions 

The MSPs in 47 Tuc are primarily very soft x-ray sources, consistent with 
thermal emission from the pulsar polar caps. The MSPs in both 47 Tuc and 
NGC 6397 seem to have a less efficient conversion of rotational spindown energy 
(P) into soft x-rays (La,) than most field MSPs, even those with correspondingly 
low values of their magnetic field at the light cylinder. For hx (X E the 47 Tuc- 
NGC 6397 samples are fit by /? = 0.5 ±0.15 whereas the updated (GCHOl) field 
(and M28) sample are consistent with the value /3 ~ 1 — 1.4 found previously 
(BT97; PCC and references therein). 

The soft colors for the 47 Tuc MSPs and the extended emission from 6397- 
A (GCHOl) indicate a lack of beamed magnetospheric emission, which may 
suggest they have different Pic values or configurations. We speculate their B 
fields may have a multipole field geometry (and thus lower dipole component) 
and be less efficient particle accelerators. Unlike MSPs in the field, those in 
dense cluster cores have a possibility of being driven back into contact and an 
accretion phase, as a re-cycled LMXB (from a MSP)! Renewed accretion (and 
6397-A is presently close to filling its Roche lobe; cf. DPM) would likely continue 
the B-field burial process thought to be responsible (e.g. Romani 1990) for field 
decay from the 10^^ G fields at NS birth to the ^ 10® G values typical of MSPs. 
The 47 Tuc MSPs have spent their r 1 Gyr lifetime in a dense (n ~ 10^ pc~^) 
cluster core, where they undergo scattering interactions with both single stars 
and other binaries. Such scattering causes the binding energy x of a hard binary 
to secularly increase. Thus, while angular momentum transfer to the secondary 
acts to detach the binary, scattering events tend to drive it back toward contact. 

A complication in this picture is that since MSPs are extremely hard bina¬ 
ries, the secular increase in x is largely the result of infrequent strong scattering 
events, which are likely to eject the binary from the cluster core or even the 
cluster. Nevertheless, it appears plausible that a typical old (r 1 Gyr) MSP 
in a dense cluster core might undergo one or more MSP-LMXB transformation 
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cycles. Thus a few percent of the MSPs at any one time may be in (or near) 
this recurrent LMXB phase. This might explain the puzzling luminous qLMXBs 
X5 and X7 (GHEOla and Heinke et al. 2002) in 47 Tuc: they might be recently 
(~ 10^“® yr) detached LMXBs, in which their underlying MSP nature is hidden 
by their evaporating (from MSP-driven winds?) fossil disks, thus explaining 
their lack of any detectable low-level accretion signatures. 

The younger system in NGC 6397 has the advantages of having been re¬ 
cently scattered out of a still higher density (~10® pc“^) core collapsed clus¬ 
ter core and possibly having exchanged its companion (Ferraro et al. 2001, 
GHEOlb). Either or both would likely have restored an accretion phase. Thus 
6397-A need not be just “born,” as suggested by Ferraro et al. (2001); it may 
instead have just been reborn. In contrast, the MSP in M28 is both ~10x 
younger, single and in a lower density core and so is unlikely to have gone 
through a renewed accretion phase. 
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